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ABSTRACT 


An  analysis  is  presented  of  the  full  scale,  free  flight  dynamic  behavior  of  a 
tri-partite  bomb  of  split-skirt  and  variable  drag  configuration. The  angular 
data  was  obtained  from  internal  instrumentation  from  the  drop  at  the  Woomera 
Test  Range  in  Salisbury,  Australia  by  the  Australian  Weapons  Research 
Establishment.  This  test  was  part  of  a  joint  research  program  on  instrumented 
bombs  undertaken  by  the  United  Kingdom,  Australia, Canada,  and  the  United 
States.  The  analysis  was  carried  out  by  the  Department  of  Aerospace  Engineer¬ 
ing,  University  of  Notre  Dame.  The  bomb  was  observed  to  experience 
Catastrophic  Yaw  in  the  early  phase  of  the  flight,  while  a  Magnus  instability 
was  evident  in  the  latter  phase.  An  analysis  of  the  last  half  of  the  flight  data 
yielded  an  excellent  determination  of  the  aerodynamic  stability  coefficients: 
cMa’CMq +CM6>  cMpcr  •  An  effort  was  also  made  to  investigate  the  motion 
of  the  bomb  as  it  passtfa  through  resonance.  Roll  lock-in  was  observed  at  this 
time  coupled  with  a  substantial  increase  in  die  magnitude  of  the  complex  angle 
of  attack,  thus  indicating  the  possibility  of  induced  side  moment  effects  and 
Catastrophic  Yaw.  Fits  of  the  Unified  Linear  Aeroballistic  Theory  to  the  data 
were  successfully  accomplished  for  this  phase  of  the  flight.  An  analysis  of  the 
motion  yielded  stability  parameters  which  indicated  the  presence  of  an  induced 
side  moment  at  resonance,  in  this  phase  of  this  analysis,  however, die  numeri¬ 
cal  procedure  used  in  fitting  the  data  caused  concern  in  that,  at  resonance, 
neither  the  Nutation  Vector  nor  the  Trim  Vector  are  rotating  (i.e.  in  body 
axes). 


This  document  is  subject  to  special  export  controls  and  each  transmittal 
to  foreign  governments  or  foreign  nationals  may  be  made  only  with  prioi 
approval  of  the  Air  Force  Armament  Laboratory  (ATBR),  Eglin  AFB, 
Florida  32S42. 
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SECTION  I 


INTRODUCTION 


At  the  present  time,  it  is  well  known  that  cruciform  finned  missiles 
and  bombs  demonstrate  three  basic  types  of  instabilities:  Magnus ,  Resonance, 
and  Catastrophic  Yaw.  During  the  past  few  years  at  the  Aero-Space 
Engineering  Department  of  the  University  of  Notre  Dame,  extensive  research 
has  been  carried  out  on  the  dynamic  behavior  of  cruciform  finned  missiles 
in  order  to  develop  a  more  fundamental  understanding  of  these  flight 
phenomena.  This  program  has  included  both  wind  tunnel  experimentation1 
and  free  flight  data  analysis.1 

Late  in  1960,  a  joint  research  program  on  instrumented  bombs  was 
formulated.  Engaged  in  this  study  are  the  British  Royal  Aircraft  Establish¬ 
ment  (RAE)  and  the  Australian  Weapons  Research  Establishment  (WRE).  At 
an  early  stage  in  the  experimental  program,  contribution  was  made  by  the 
United  States  Navy,  Air  Force,  and  Army  in  the  form  of  wind  tunnel  and 
computer  facilities P  In  this  joint  program,  effort  was  directed  toward 
investigation  of  the  dynamic  characteristics  of  less  conventional  stabilizers. 
In  particular,  a  split-skirt  tail  configuration  was  used  for  experimentation 
in  the  hope  that  it  would  eliminate  one  or  more  of  the  instabilities  character¬ 
istic  to  cruciform  finned  missiles.  However,  free  flight  tests  of  this  config¬ 
uration,  performed  at  the  Woomera  Test  Range  in  Salisbury,  Australia4, 
indicated  the  possibility  of  these  instabilities. 

It  was  felt  that  application  of  the  techniques  developed  at  the  University 
of  Notre  Dame  would  br  useful  in  examining  the  flight  characteristics  of  this 
new  bomb  configuration.  Thus,  the  purpose  of  this  report  is  to  analyze  the 
dynamic  behavior  of  one  test  round  of  a  split-skirt  bomb.  Approximately 
1000  points  of  telemetry  records  of  the  angular  orientation  data  as  a  function 
of  time  were  fitted  to  the  Unified  Linear  Aeroballistic  Theory  by  the  Method 
of  Differential  Corrections  as  described  in  Reference  2.  From  the  results 
of  this  fitting  technique,  the  following  aerodynamic  stability  coefficient 
derivatives  were  determined:  pitching  moment,  Cm*  ;  pitch  damping 
moment, C*u  +Cm*  ;and  Magnus  moment,  Cmp«  •  In  addition,  an 
extension  of  me  Unified  Linear  Aeroballistic  Theory  was  made  in  order  to 
examine  the  effects  of  me  induced  side  moment  coefficient  derivative, 

Ciuy ip. fciN+b* ,  at  resonance.  The  results  of  mis  portion  of  the  analysis, 
however,  are  questionable  due  to  a  combination  of  die  flight  characteristics 
at  resonance  and  die  numerical  procedures  employed  in  fitting  the  data. 
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SECTION  II 


THEORY 


For  the  sake  of  completeness,  a  development  of  the  Six  Degree  of 
Freedom  Equations  of  Motion,  as  well  as  the  Unified  Linear  Aeroballistic 
Theory  and  its  non-linear  extension  will  be  presented. 


Six  Degree  of  Freedom  Equations  of  Motion 


Consider  two  sets  of  orthogonal  axes:  a)  space -fixed  axes  x,y,z;  and 
b)  missile-fixed  axes  X,Y,Z.  Beginning  with  Newton’s  Laws  of  Translational 
and  Angular  Motion,  it  is  proposed  to  obtain  a  set  of  equations  of  motion  for 
a  missile  with  six  degrees  of  freedom,  whose  solution  will  yield  the  criteria 
for  dynamic  stability  with  expressions  for  the  pertinent  aerodynamic  forces 
and  moments.  It  should  be  noted,  however,  that  the  aerodynamic  parameters 
are  with  respect  to  the  set  of  axes  fixed  to  the  missile  and  that  Newton's 
Laws  of  motion  are  valid  only  for  a  set  of  axes  which  are  stationary  or 
translate  with  a  constant  velocity.  According  to  Goldstein*,  however,  the 
following  transformation  from  the  missile-fixed  axes  to  the  space-fixed  axes 
may  be  applied: 


space  fixed 

Where  £1  is  the  angular  velocity  of  the  coordinates  of  the  wobbling  axes. 


si 


wobbling  axes 


+XiXV 


(1) 


Consider  first,  Newton's  Law  of  Translational  Motion  applied  to  the 
space -fixed  axes. 

t-  (2) 


where 


The  components  of  \/  are  the  linear  velocity  components  with  respect  to 
the  missile-fixed  axes.  Noting  that  mass  is  constant,  application  of  Eq.(l) 
toEq.  (2)  yields 


where 


2 


Jfl= 


-Ax 

Ay 

L^a. 


,  the  angular  velocity  components 
of  the  coordinate  system. 


Thus, 


(3) 


In  Aeroballistic  axes* 

2t 


where 


rV7vT 


(4) 


Thus,  performing  the  indicated  differentiation  and  vector  operation  in  Eq. 
(3)  and  writing  the  force  in  terms  of  its  components , 


FX(  =  X),  FY  (  -  Y),  FZ  (  =  Z)  yields 

X  a  jtor-irnr) 

Y  =  Tmar+miru 

2s>nur-Y«^u 

Consider  Newton's  Law  of  Angular  Motion  applied  to  the  space-fixed 


(5) 


axes, 


(6) 


where 


x= 


r  £«»  Xx« 
Xrx  Xyy  Xyi 

1**  I*y  X»I 


the  moment  of  inertia  tensor, 


•Aeroballistic  axes  are  orthogonal  axes  fixed  to  the  missile  but  which  do  not 
roll  with  it.  The  origin  is  fixed  at  the  center  of  gravity  and  the  X  axis  is 
coincidental  with  the  symmetric  axis  of  the  missile. 
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and, 


CO  = 


f 

% 

r. 


the  angular  velocity  components 
in  the  missile -fixed  axes. 


Assuming  the  principal  axis  is  coincident  with  the  geometric  axis,  the 
moment  of  inertia  tensor  may  be  diagonalized,  and  assuming  rotational  mass 
symmetry,  Ir=  Xgs  X  •  Thus, 

i*  o  6 
1=  o  i  o 
o  o  : 

Application  of  Eq.  (1)  to  Eq.  (6)  yields 

or  expanding  into  components  for  aeroballistic  axes 


’I*  o  off* 

■  m 

o 

'lx 

O 

mm 

o 

OX 

+ 

5 

O 

I 

o 

Lo  o  xJLrJ 

Jr 

O 

o 

I 

t. 

(7) 


Performing  the  indicated  vector  operations  in  Eq.  (7)  and  writing  the  moment 
in  terms  of  its  components  Mx(  *  L),  My(  =  M),  Mz(  *  N)  yields 


M=ii+w*r 


(8) 


Thus,  Eqs.  (5)  and  (8)  represent  the  differential  equations  of  motion  In 
aeroballistic  axes  for  a  free  flight  missile  having  six  degrees  of  freedom. 


Unified  Linear  Aeroballistic  Theory 

The  solution^  of  these  equations  for  the  complex  angle  of  attack,  subject 
to  die  constraints  of  constant  roll  rate  and  velocity,  is  the  following: 


(9) 


where 


<*A  =  complex  angle  of  attack  (aeroballistic) 
T^=  nutation  vector  (aeroballistic) 

=  precession  vector  (aeroballistic) 

T<Ta  =  trim  vector  (aeroballistic) 

yaw  of  repose  vector  (aeroballistic) 


(10) 


^=#d±4) 


an 


T=-J 


\FT 


(12) 


c- 

4IAdQCM. 


(13) 


The  solution  for  the  complex  angle  of  attack  may  be  expressed  in  body 
axes*  by  application  of  the  transformation 


S.scxxe 

*  Body  axes  are  aeroballistic  axes  which  roll  with  the  missile. 


(14) 
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Thus,  application  of  Eq„  (14)  to  Eq.  (9)  yields  the  solution  in  terms  of  body 
axes: 


where 


complex  angle  of  attack  (body) 
nutation  vector  (body) 

precession  vector  (body) 

trim  vector  (body) 

yaw  of  repose  vector  (body) 


Recognizing  the  fact,  however,  that 
systems  is  the  rolling  velocity,  or 


die  onljkdifferepce  between  the  two  axis 

*■  j  «te. 


CONb=  60^-  Jp  (16) 

and, 

<17> 

Thus,  die  solution  for  the  complex  angle  of  attack  in  body  axes  may  be 
written  as 

where  die  vector  dynamic  damping  factors  and  frequencies  are  defined  by 
Eqs.  (10)  through  (13),  (16)  and  (17). 

The  necessary  and  sufficient  conditions  for  the  dynamic  stability  of  a 
fin  stabilized  missile  are  as  follows: 


1.  The  gyroscopic  stability  factor  must  be  negative,  s<0,  and, 

2.  The  nutation  and  precession  dynamic  damping  factors  must 
be  negative,  'X^p<  O  • 
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Catastrophic  Yaw  Theory^ 

A  unique  instability,  not  accounted  for  by  the  Unified  Linear  Aero- 
ballistic  Theory^,  was  occasionally  observed  in  the  trajectories  of  a  certain 
fin- stabilized  missile.  This  instability  was  characterized  by  two  phenomena: 
1)  failure  of  the  missile  to  attain  its  steady  state  rolling  velocity,  and  2)  a 
rapid  increase  of  the  pitching  and  yawing  motion  to  extremely  large  angles  of 
attack.  These  phenomena  are  labeled  "Roll  Lock-In"  and  "Catastrophic  Yaw," 
respectively.  7,8,9 

Initially,  the  rolling  velocity  increases  due  to  fin  cant,  but  when  it 
attains  a  value  equal  to  the  aeroballistic  nutation  frequency,  it  "locks  In"  at 
this  particular  value  rather  than  continuing  to  seek  its  steady  state  rate, 
which  is  much  larger.  With  respect  to  the  pitching  and  yawing  motion,  how¬ 
ever,  initially  the  magnitude  decreases;  but  when  the  rolling  velocity  locks 
in,  the  motion  increases  to  extreme  values  which  may  even  cause  the  missile 
to  tumble. 

When  observed  in  special  wind  tunnel  and  full  scale  free  flight  tests, 
the  phenomenon  of  roll  lock-in  appears  to  be  traceable  to  the  influence  of 
the  non-linear  induced  roll  moment.  11  Indeed,  when  this  moment  is 
introduced  into  the  classical  theory  for  pure  rolling  motion,  a  method  for 
predicting  the  critical  angle  of  attack  and  roll  trim  angle  is  available.  Of 
major  concern  here,  however,  is  the  Catastrophic  Yaw. 

As  established  by  Nicola  ides®  and  supplemented  by  flight  performance 
data,  the  catastrophic  growth  of  the  pitching  and  yawing  motion  occurs  when 
the  missile  is  in  "lunar  motion,"  that  is,  when  the  angle  between  a  reference 
fin  and  the  plane  of  complex  angle  of  attack  is  constant.  The  Unified  Linear 
Aeroballistic  Theory  predicts,  for  a  statically  stable  missile,  that  this  type 
of  motion  may  exist  for  any  of  three  cases:  1)  pure  nutational  motion  in 
which  the  roll  rate  is  equal  to  the  aeroballistic  nutation  frequency,  2)  pure 
trim  in  which  lunar  motion  exists  for  all  values  of  the  rolling  velocity,  and 
3)  a  combination  of  bod)  in  which  die  rolling  velocity  is  equal  to  the  aero¬ 
ballistic  nutation  frequency.  The  Linear  Theory,  applied  to  pure  nutational 
motion,  fails  to  account  for  die  catastrophic  growth  observed.  Also, 
maximum  amplification  of  the  trim  occurs  when  the  rolling  velocity  equals 
the  aeroballistic  nutation  frequency,  which  is  known  as  "Resonance  Instability. 
Initially,  one  would  suspect  this  phenomenon  to  be  the  cause  of  die  observed 
motion;  this  amplification,  however,  may  be  extracted  from  the  linear 
theory,  which  also  yields  approximate  values  when  non-linearities  exist  in 
the  force  and  moment  system.  In  neither  case  does  die  linear  theory  account 
for  die  observed  catastrophic  growth  of  complex  angle  of  attack.  Thus,  it  is 
evident  that  a  further  understanding  and  evaluation  of  die  fluid  forces  and 
moments  acting  on  die  missile  is  necessary. 


While  the  dependence  on  roll  orientation*  is  automatically  d«tect«d 
from  the  linear  theory,  its  essential  role  in  the  explanation  of  roll  lock-in 
suggests  that  there  may  be  a  force  and  moment  associated  with  it.  Indeed, 
wind  tunnel  tests  do  reveal  two  additional  effects  of  roll  orientation:  1)  the 
normal  force  and  its  moment  are  modified,  and  more  important,  2)  a  side  > 

force  and  moment  are  found  to  exist.  The  side  moment  was  observed  to  vary 
sinusoidally  with  the  roll  orientation. 

Thus,  in  addition  to  the  forces  and  moments  involved  in  the  linear 
theory,  the  following  are  now  introduced  (for  a  cruciform  finned  missile): 


SIN 4*  S  A  SIN4* 

and 


(19a) 


M(y4pO=  My4<x(X  SIN^tfsCM^P^QAd  SIN4Jf  (19b) 

In  order  to  determine  the  contribution  of  these  forces  and  moments  to  the 
dynamic  stability  of  a  missile,  they  must  be  added  to  the  classical  aero¬ 
dynamic  system.  A  general  solution  for  the  complex  motion  involving  these 
new  terms  is  not  possible;  however,  according  to  Nicolaides^,  it'  the 
characteristic  motion  is  assumed  to  be  "lunar”  an  approximate  solution  is 
possible.  Employment  of  the  perturbation  approach  to  this  type  of  motion, 
as  presented  in  Reference  6,  yields  a  variation,  due  to  the  additional  forces 
and  moments,  in  the  equation  for  the  dynamic  damping  factors.  Thus, 
neglecting  changes  in  the  normal  force  due  to  roll  orientation,  this  equation 
assumes  the  following  form: 

*  (20) 

Comparison  with  the  resulting  dynamic  damping  factors  from  linear  theory, 
Eq.  (10),  ami  denoting  them  as  ,  yield  the  following  relation: 

(SIN4¥)X  (21) 


♦Roll  orientation  angle,  )f  ,  is  defined  as  the  angle  between  a  reference  fin 
and  the  plane  formed  by  the  complex  angle  of  attack. 
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It  should  be  noted  that  while  in  stability  derivative  notation,  a  subscript 
implies  partial  differentiation  with  respect  to  the  variables,  is  not 

a  derivative  with  respect  to  y  but  a  function  of  it.  Thus,  it  ia  completely 
consistent  with  the  notation .Cm^  — ♦  C ( y)  .  Hence,  in  Catastrophic 
Yaw  Theory  the  necessary  and  sufficient  conditions  for  the  dynamic  stability 
of  a  fin^stabilized  missile  are  as  follows: 

1.  The  gyroscopic  stability  factor  must  be  negative,  s  <  0,  and, 

2.  The  nutation  and  precession  vector  dynamic  damping  factors  must 

be  negative,  O  • 
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SECTION  III 


COMPUTATIONAL  PROCEDURE 


The  computer  program  of  Reference  12  was  used  to  fit  the  complex 
angular  data4  to  the  quadricyclic  equation,  Eq.  (18),  by  the  Method  of 

Differential  Corrections 1 3, 14 #  The  stability  parameters ,  .CJWpp*.  > 

»  were  determined  as  a  function  of  time  from  release  by  fitting 
the  Unified  Linear  Aeroballistic  Theory  to  small  segments  of  the  data  and 
overlapping  these  segments. 


Two  assumptions  are  made  in  using  this  technique  to  obtain  the 
stability  parameters. 

1.  The  segment  of  fit  is  chosen  small  enough  to  insure  that  the 
total  velocity  and  roll  rate  are  essentially  constant  in  the 
time  interval  defining  the  segment. 

2.  Enough  points  of  the  complex  angular  data  are  fitted  in  each 
of  these  segments  to  well  define  the  stability  parameters  for 
the  time  interval  defining  that  segment. 

The  fact  that  these  segments  are  overlapped  enables  one  to  obtain  the 
stability  parameters  as  continuous  functions  of  time  from  release. 


If  the  complex  angular  motion  is  such  that  both  the  Nutation  and 
Precession  Vectors,  their  damping  factors  and  frequencies  may  be  obtained 
from  the  fitting  procedure;  the  aerodynamic  stability  derivatives,  Cm*  * 
Cm*  4 Cm*  >  and  Cnf  may  be  computed  as  found  in  Reference  2. 

This  general  procedure  is  applied  in  the  following  manner: 


1.  Substituting  Eqs.  (16),  (17),  (18),(13)#(12)  into  Eq.  (11)  and 
solving  for  C**  yields, 

CKVcii^7l[  «^.**W*V  <»> 

2.  Using  predetermined  values  for  Cg*  »  Cm^+Cm*  may  be 

computed  by  adding  the  expressions  for  and  \f  ,  and 
solving  for  yields, 

<VC*«4Srt.-C"«  +3&<Xh4Xp>]  (23) 
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3.  Having  solved  for  *  solving  either  the  X*  or  Xp 

expression  for  yields , 

(24) 

(Cm,*Cm^Ki it)J 

If,however,  the  complex  angular  motion  is  such  that  the  Precession 
Vector  is  zero  or  is  small  that  the  fitting  routine  is  not  able  to  determine 
it  within  a  reasonable  degree  of  accuracy, 

1.  Substituting  Eqs.  (16),  (13),  and  (12)  into  the  expression  for 
Eq.  (11),  and  solving  for  Cm.,  yields, 

\ 

'*  (25> 

Since  in  this  case,  Xp  is  not  known,  an  analytical  solution  for 
Ck  is  not  possible.  Thus;  in  order  to  calculate  CWp*  » 

Ck.  must  be  known.  If  this  is  the  case,  Cm,*  may 

stilTbe  calculated  as  follows:  r 

2.  Solve  the  Xm  expression  for  Cm^  yielding 

**»~c*«  0-x)  -f^«b,«CO0«^26) 

In  the  Catastrophic  Yaw  Theory,  consider  Eq.  (21), 

=)!*,  if^CM^tSIM+VJX  (21) 

solving  for  the  induced  side  moment  term, 


3&c"**ts,M4W'c“  x"'x" 


(27.) 


(27b) 


n 


and  making  use  of  Eq.  (19b), 


(28a) 

(28b) 

may  be  obtained 

(29a) 


(29b) 


SECTION  IV 


PRELIMINARY  ANALYSIS 

Configurational  Description  and  Free  Flight  Conditions 


The  configuration  analyzed  is  a  free  flight,  variable  drag  bomb  with  a 
split-skirt  tail,  which  has  a  skirt  opening  of  10  degrees.  A  schematic  is  pre¬ 
sented  in  Fig.  1.  The  mass  parameters  for  this  test,  Round  739,  are  pre¬ 
sented  in  Table  1 . 

This  bomb  was  dropped  from  an  altitude  of  45000  feet;  the  time  of  fall 
was  recorded  as  67.47  seconds.  Angular  orientation  data  was  determined 
from  .002  second  to  63.  39  seconds  from  release  at  intervals  of  .047 
second. 

The  angular  orientation  data4  for  this  analysis  consists  of  the  two 
components  of  the  complex  angle  of  attack  \/8+Loc  )  in  body  axes .  These 
components  are  presented  in  Figs.  2  and  3  as  a  function  of  time  from  re¬ 
lease,  and  the  magnitude  of  the  complex  angle  of  attack,  l#  l  ,  is  shown 
in  Fig.  4,  while  the  complex  motion  is  presented  in  Fig.  5.  This  report 
consists  of  a  dynamic  stability  analysis  from  12.05  to  62. 13  seconds  from 
release. 

The  following  trajectory  parameters  for  Round  739  are  presented  in 
Fig.  6  as  a  function  of  time  from  release:  roll  rate,  velocity,  Mach  number, 
dynamic  pressure,  and  Reynolds  number  per  foot.  Upon  preliminary  exami¬ 
nation  of  this  data,  it  was  noted  that  the  roll  rate,  while  increasing  initially, 
levels  off  at  approximately  12  seconds  until  20  seconds  at  which  time  it 
speeds  up  towards  its  steady  state  value  at  the  end  of  the  flight.  This  phenom¬ 
enon,  combined  with  the  observed  magnitude  of  complex  angle  of  attack  At 
approximately  the  same  time,  indicated  the  possibility  of  Roll  Lock-In  and 
Catastrophic  Yaw  effects  in  the  early  phase  of  the  flight.  It  was  further 
noted  that,  although  Mach  number  variations  were  not  extremely  large,  the 
flight  was  characterized  by  transonic  velocities  from  approximately  35 
seconds  to  42  seconds  from  release. 

From  the  magnitude  of  the  complex  motion,  two  dominant  character¬ 
istics  of  the  motion  were  observed:  a)  the  rapid  increase  and  decrease  of 
the  magnitude  of  complex  angle  of  attack  at  approximately  14  seconds  from 
release,  and  b)  the  characteristically  undamped  motion  from  25  seconds  from 
retease  to  die  end  of  the  flight.  Referring  these  observations  to  the  complex 
motion,  Fig.  5,  it  was  noted  that,  at  approximately  14  seconds  from  release, 
one  of  the  vector  frequencies  changed  sign,  i.e. ,  one  of  the  "loops," 
characteristic  to  the  motion,  reversed  direction.  Since  the  motion  is  in  body 
axes,  a  change  in  sign  of  the  body  nutation  frequency  would  indicate  that  the 
bomb  had  passed  through  resonance  at  this  time.  Also,  from  Fig.  5  it  was 
noted  that  from  25  seconds  from  release  to  the  end  of  the  ilight,  the  motion 
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TABLE  I.  MASS  PARAMETERS  OF  SPLIT-SKIRT  BOMB 


Release  Weight  (lb) 

939.0 

Length  (in.) 

144.0 

Maximum  Body  Diameter  (in. ) 

18.75 

Axial  Moment  of  Inertia  (slug-ft^) 

8.8 

Transverse  Moment  of  Inertia  (slug-ft^) 

171.3 

C.  G.  Position  from  Nose  (%  body  length) 

30.7 
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is  characterized  by  undamping  corcentric  circular  motion  indicative  of  a 
Magnus  instability.  Thus,  preliminary  investigation  of  the  angular  orienta¬ 
tion  data  indicated  that  two  of  the  three  basic  types  of  instabilities  character¬ 
istic  to  cruciform  finned  missives  might  be  present:  Resonance  and  Magnus 
Instability. 


Fitting  Procedure 

The  fitting  technique  of  Reference  12  is  organized  in  such  a  manner  as 
to  allow  considerable  flexibility  in  its  application  to  various  sets  of  data.  In 
the  case  of  six  degree  of  freedom  data,  the  motion's  characteristics,  such  as 
rate  of  change  of  the  magnitude  of  the  angle  of  attack  and  changes  in  the  shape 
of  the  complex  motion  (circular,  elliptical,  etc.)  generally  dictate  the  pro¬ 
cedures  to  be  followed  to  obtain  the  best  possible  fit.  Another  important 
factor  which  must  be  considered  is  the  possibility  of  error  in  the  data  and 
how  this  error  may  affect  the  results.  Both  factors  were  investigated  in  great 
detail  prior  to  application  of  the  fitting  routine. 

In  order  to  detect  possible  error  in  the  data,  investigation  was  directed 
towards  the  magnitude  of  the  complex  motion,  Fig.  4,  in  which  considerable 
oscillations  were  observed.  Up  to  approximately  25  seconds  from  release, 
at  which  time  the  complex  motion  assumes  a  circular  pattern,  it  was 
assumed  that  these  oscillations,  which  decreased  in  amplitude  with  time, 
were  due  to  the  superposition  of  a  damping  vector  with  the  vector  dominant 
in  the  motion.  After  25  seconds  from  release,  however,  the  motion  is 
circular  about  a  center  at  nearly  zero  complex  angle  of  attack,  indicating 
that  the  motion  is  probailv  monocyclic  with  a  small  trim.  Thus,  it  was 
thought  that  the  observed  oscillations  might  be  due  to  inherent  error  in  the 
data.  Rhodes  and  Shannon^  indicate  that  the  differential  pressure  incidence 
meter  used  to  measure  the  angular  orientation  data  has  phase  lags  associated 
with  the  tube  lengths  leading  from  the  surface  tapping  points  to  the  pressure 
transducers.  These  phase  lags  are  estimated  in  Reference  15.  From 
previous  experience  with  the  fitting  routine  ,  it  was  thought  that  these  oscil¬ 
lations  might  influence  the  determination  of  the  dynamic  damping  factors. 

The  procedure  to  be  followed  in  fitting  the  angular  orientation  data  was 
heavily  dependent  on  the  general  characteristics  of  the  complex  motion. 

From  12.05  seconds  from  release  to  approximately  16.6  seconds,  the  com¬ 
plex  motion  consists  of  large  "loops."  At  this  time,  the  "loops"  become 
much  smaller  and  begin  to  damp  out  while  the  overall  motion  is  elliptical. 

At  24.5  secondi  the  motion  is  of  a  circular  nature  and  continues  in  this 
manner  until  the  end  of  the  flight.  The  angular  orientation  data,  tabulated 
at  every  0.047  second,  was  fitted  to  the  quadricytlic  equation,  Eq.  (18), 
over  various  numbers  of  cycles  of  the  complex  motion,  depending  on  what 
section  of  the  data  was  being  fitted.  Beginning  at  16.62  seconds,  the  fit  was 
receded  in  time  in  increments  of  0.047  second,  each  sectional  fit  encom- 


passingO.  Scyclc ,  until  the  effects  of  the  release  disturbance  were  felt  at 
12.05  seconds  and  the  fit  was  terminated.  The  fit  was  then  advanced  from 
16.62  seconds  at  increments  of  0.  13  second,  each  sectional  fit  encompassing 
0.75  cycle,  until  24.5  seconds  at  which  time  the  Precession  Vector  had 
damped  to  such  an  extent  that  the  fitting  routine  could  not  determine  it,  its 
dynamic  damping  factor,  and  frequency  to  within  a  reasonable  degree  of 
accuracy.  In  the  previous  two  sections  of  data  fitted,  the  character  of  the 
motion  was  changing  rapidly,  thus  dictating  that  each  fit  encompass  only  a 
fraction  of  a  cycle.  In  this  case  the  number  of  points  per  cycle  was 
sufficiently  large  to  insure  a  good  fit. 

At  24.5  seconds,  initial  conditions  with  sufficient  accuracy  to  commence 
the  fit  were  not  able  to  be  extracted  from  the  data.  At  34.97  seconds, 
however,  initial  conditions  were  obtained,  and  the  fit  was  receded  in  time  to 
24.5  seconds  at  increments  of  0. 141  second,  each  sectional  fit  encompassing 
2  cycles.  The  ability  to  obtain  an  accurate  fit  with  a  larger  number  of  cycles 
is  due  to  the  fact  that  the  character  of  the  motion  was  not  changing  too 
severely  from  24.5  seconds  to  34.97  seconds.  From  34.97  seconds  to  the 
end  of  the  flight,  the  motion  is  essentially  circular.  Thus,  more  cycles  may 
be  included  per  sectional  fit  without  loss  of  accuracy.  This  increase  in  the 
number  of  cycles  is  also  necessary  since  the  number  of  points  per  cycle  was 
decreasing.  Therefore,  from  34.97  seconds  from  release  to  the  end  of  the 
flight,  the  fit  was  advanced  0.23  second,  each  sectional  fit  encompassing  3 
cycles,  until  the  fit  was  terminated  at  62. 13  seconds  from  release.  A 
summary  of  the  parameters  and  techniques  used  in  the  fitting  procedure  is 
presented  in  Table  II. 

Thus,  by  overlapping  these  sectional  fits  the  stability  parameters  of 
Eq.  (18)  were  determined  as  a  function  of  time  from  release.  Using  position 
coordinate  data  and  the  magnitude  of  the  complex  angle  of  attack,  the 
functional  dependence  of  the  stability  derivatives  with  Mach  number,  Reynolds 
number,  and  angle  of  attack  was  analyzed. 


TABLE  II.  FITTING  ROUTINE  AND  PARAMETERS 


Initial 
me  (sec) 

Final 

Time  (sec) 

Time 

Increment  (sec) 

Number  of  Cycles 
per  Sectional  Fit 

16.62 

12.05 

-.047 

.5 

16.62 

24.5 

.  13 

.75 

34.97 

24.5 

-.141 

2.  0 

34.97 

62. 13 

.23 

3.0 

37 
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SECT  ION  V 


ANALYSIS  OF  RESULTS 


The  fit  of  the  complex  angular  motion  to  the  quadricyclic  equation, 

Eq.  (18),  yielded  a  total  probable  error  which  ranges  from  0. 12  degree  to 
2.3  degrees.  Fig.  7  shows  the  probable  error  of  fit  (Figs.  7a  to  7c)  and  the 
percent  error  of  fit  (Figs.  7d  to  7f)  as  a  function  of  time  from  release .  The 
percent  error  is  with  respect  to  the  magnitude  of  the  vectors  fitted.  It  is 
seen  that  the  maximum  percent  error  obtained  is  approximately  1%  of  the 
total  motion,  which  indicates  a  good  fit  of  the  theory  to  the  data.  At  24.5 
seconds  from  release,  there  is  a  sharp  increase  in  the  probable  erroi  from 
0.2  degree  to  0.6  degree.  There  are  two  reasons  for  this  discontinuity: 

1)  At  the  end  of  the  previous  fit  (24.5  seconds)  the  magnitude  of  the 
Precession  Vector  had  damped  to  such  an  extent  that  the  fitting  routine  was 
no  longer  able  to  determine  it,  its  dynamic  damping  factor  and  frequency. 
Hence,  at  the  beginning  of  the  next  fit,  these  small  values  of  the  Precession 
Vector  appeared  in  the  residuals  as  error  until  the  vector  damped  to  a 
negligible  amplitude  at  approximately  26.5  seconds  from  release.  2)  A  higher 
probable  error  at  the  beginning  of  a  fit  has  been  found  to  be  characteristic 
to  the  fitting  procedure  due  to  the  fact  that  it  takes  several  iterations  for  the 
initial  approximations  to  attain  sufficient  accuracy.  Thus,  a  combination  of 
these  observations,  one  characteristic  to  the  data,  the  other  to  the  fitting 
procedure,  led  to  the  discontinuity  observed  in  the  probable  error.  The 
amplitude  of  the  probable  error,  indicating  an  extremely  good  fit  to  the 
complex  angular  data,  implied  generally  accurate  fits  to  the  linear  stability 
parameters  of  the  quadricyclic  equation. 


Linear  Stability  Parameters  and  Probable  Errors 

In  fitting  the  complex  angular  motion,  Fig.  5,  to  the  quadricyclic 
equation,  Eq.  (18),  it  was  found  that  the  Nutation  Vector  and  Trim  Vector 
were  present  throughout  the  entire  flight.  The  Precession  Vector,  however, 
was  not  fitted  after  24.5  seconds  from  release,  while  the  Yaw  of  Repose 
Vector  was  found  to  be  negligibly  small  throughout  the  entire  flight. 


Nutation  Vector  and  Frequency 

Figs.  8  and  9  show  the  magnitude  of  Nutavion  Vector  and  its  probable 
error  as  a  function  of  time  from  release.  The  fit  yielded  an  accuracy  which 
ranged  from  0.04  degree  to  l.l  degrees.  It  should  be  noted  that  the  largest 
probable  error,  1. 1  degrees,  occurred  when  the  magnitude  was  approxi¬ 
mately  33  degrees.  Comparison  of  the  magnitude  of  k*  and  the  magnitude 
of  the  complex  motion  from  25  seconds  until  the  end  of  the  flight  indicates  that 
practically  the  entire  motion  was  pure  nutation.  The  slight  oscillations 
appearing  in  the  magnitude  are  felt  to  be  due  to  the  inherent  error  in  the 
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data,  as  discussed  in  the  Preliminary  Analysis. 

The  nutation  frequency,  in  body  axes,  and  its  probable  error  are  pre¬ 
sented  &s  a  function  of  time  from  release  in  Figs.  10  and  11.  Fig.  10  shows 
an  approximately  linear  increase  in  from  0.6  rad/sec  at  14.7 

seconds  to  16.  3  rad/sec  at  57  seconds,  at  which  time  it  remains  constant 
until  61  seconds  when  a  slight  decrease  is  observed.  The  results  of  Fig.  11 
indicate  an  extremely  accurate  fit.  This  was  expected,  however,  from 
previous  experience  with  the  fitting  procedure  which  indicated  that,  in 
general,  the  vector  frequencies  were  the  most  accurately  determined  para¬ 
meters.  Fig.  11  shows  a  range  of  probable  errors  from  .007  rad /sec  to 
0.075  rad /sec.  In  this  case,  however,  the  maximum  probable  error  did  not 
occur  when  the  frequency  was  the  largest.  On  the  contrary,  the  frequency 
was  changing  sign  at  this  time.  From  Eq.  (16)  it  is  obvious  that  at  this  time 
the  aeroballistic  nutation  frequency  is  equal  to  the  rolling  velocity,  i.e.,the 
bomb  is  at  resonance.  This  is  in  complete  agreement  with  the  observations 
made  in  the  Preliminary  Analysis. 

The  largest  probable  error  occurring  at  the  time  of  zero  body  nutation 
frequency  is  felt  to  be  of  considerable  importance.  In  the  fitting  procedure, 
the  number  of  points  to  be  fitted  per  section  is  determined  from  a  calculation 
involving  division  by  the  nutation  frequency.  Hence,  if  the  angular  data  is  in 
body  axes,  the  number  of  points  to  be  fitted  at  resonance  will  go  to  infinity. 

In  order  to  avoid  this,  a  slight  alteration  was  made  in  the  fitting  routine: 
the  number  of  points  to  be  fitted  was  held  constant  when  the  magnitude  of  the 
nutation  frequency  was  near  0  rad/sec.  This  alteration  induced  a  large 
rate  of  change  of  CO^  through  resonance  as  shown  in  Fig.  10a.  Therefore, 
in  this  small  time  range,  the  linear  theory  was  violated.  This  explains  the 
sharp  increase  in  probable  error  of  fit  obtained  at  this  time,  as  shown  In 
Fig.  7.  Similar  effects  were  noted  in  the  probable  errors  of  the  stability 
parameters  at  this  time. 


Precession  Vector  and  Frequency 

The  magnitude  of  the  Precession  Vector  and  its  probable  error  are 
presented  in  Figs.  12  and  13.  The  magnitude  ranges  from  14  degrees  at 
approximately  14  seconds  to  1  degree  at  24.5  seconds.  At  24.5  seconds  the 
vector  decreased  to  such  a  small  magnitude  that  it  was  no  longer  able  to  be 
extracted  from  the  data  as  the  Precession  Vector.  The  larger  probable  error 
of  fit  at  this  time  indicated  that  |1<^|  appeared  in  the  residuals  as  error 
until  its  magnitude  was  negligible.  The  greatest  error  in  |K^t  ,0  .62 
degree  ,  occurred  at  approximately  14  seconds  from  release,  when  the 
magnitude  of  the  vector  was  die  largest.  The  sharp  peak  in  the  probable 
error  at  this  time  is  felt  to  result  directly  from  the  violation  of  the  linear 
theory  as  discussed  previously.  The  error  is  seen  to  become  slightly  more 
random  as  the  vector  begins  to  damp  out  of  the  motion. 
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As  presented  in  Fig.  14,  the  body  precession  frequency  varies  slightly 
throughout  the  section  fitted,  ranging  from  6  rad /'sec  at  12.05  seconds  to 
11  rad /sec  at  24.5  seconds.  The  probable  error,  shown  in  Fig.  15,  increases 
from  0.03  rad/sec  at  13  seconds  to  0.25  rad/sec  at  24.5  seconds.  The  slight 
increase  observed  as  the  fit  approaches  its  termination  at  12.05  seconds  is 
felt  to  be  due  to  the  effect  of  release  disturbance.  The  increase  in  probable 
error  beginning  at  approximately  22  seconds  results  from  the  fact  that  the 
Precession  Vector  is  in  the  process  of  damping  out  at  this  time.  The  effects 
are  seen  in  the  slight  variations  of  (Op.  in  Fig.  14, 

Trim  Vector 

The  Trim  Vector  .  shown  in  Fig.  16,  is  seen  to  reach  its  maximum 
value,  21  degrees,  at  approximately  14.2  seconds,  the  resonant  condition, 
as  expected  from  the  analysis  of  the  body  nutation frequency ,  F-v  t  " '  The 
probable  error,  presented  in  Fig.  17,  indicates  an  extremely  good  fit  even 
though  the  magnitude  of  the  Trim  Vector  is  very  nearly  zero  degrees  during 
the  last  37  seconds  out  of  the  flight.  Again,  the  maximum  probable  error 
occurs  when  the  magnitude  of  the  vector  is  also  a  maximum .  As  in  the  case 
of  the  Precession  Vector,  the  sharp  increase  in  the  probable  error  occurs 
when  the  linear  theory  was  violated. 

At  approximately  24  seconds  from  release ,  when  the  complex  motion 
begins  to  assume  a  circular  nature,  the  Trim  Vector  is  almost  negligible  . 
This,  combined  with  the  magnitude  of  the  Nutation  Vector  and  the  loss  of 
the  Precession  Vector,  completely  explains  the  concentric  circular  motion 
about  an  almost  zero  degree  trim  observed  in  the  complex  plane.  A  com* 
parison  of  the  magnitude  of  the  Nutation,  Precession,  and  Trim  Vectors  at 
resonance  indicate  that  the  major  cause  of  the  large  angle  of  attack  observed 
is  the  amplification  of  the  Trim  Vector.  Fig.  18  shows  Ikt|  as  a  function 
of  die  ratio  of  the  roll  rate  to  the  aeroballistic  nutation  frequency.  The 
Trim  Vector  is  seen  to  reach  maximum  amplitude  at  the  point  where  the 
ratio  is  1. 

Nutation  and  Precession  Dynamic  Damping  Factors 

As  anticipated  in  the  Preliminary  Analysis,  the  oscillations  in  the 
angular  orientation  data  had  an  undesirable  effect  on  die  determination  of 
the  dynamic  damping  factors.  This  effect  was  predominant  when  die  magni¬ 
tude  of  the  complex  angle  of  attack  was  not  larger  than  approximately  1* 
degrees.  This  occurred  from  16. 6  seconds  to  the  end  of  die  flight  At  these 
times  the  complex  motion  was  characterized  by  quasi-circular  and  circular 
motion.  That  is,  from  16.6  seconds  to  24.5  seconds,  the  Precession  Vector 
was  in  the  process  of  damping  out,  and  from  this  time  until  the  end  of  die 
flight  the  motion  was  characteristically  nutational.  Up  to  16.6  seconds, 
however,  the  magnitude  of  the  angle  of  attack  was  either  extremely  huge, 
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or  was  changing  so  rapidly  that  the  oscillations  were  concealed.  Further¬ 
more,  the  times  at  which  the  oscillations  played  an  important  role  in  the 
determination  of  the  dynamic  damping  factors  were  characterized  by  either 
a  very  small  value  in  the  nutational  mode  or  the  termination  of  the  vector 
in  the  precessional  mode.  Thus,  since  the  dynamic  damping  factors  play 
an  important  role  in  the  determination  of  the  aerodynamic  coefficients ,  con¬ 
siderable  care  was  employed  in  evaluating  the  results. 

Fig.  19  shows  the  Nutation  Vector  dynamic  damping  factor  as  a  function 
of  time  from  release.  As  is  evident  from  this  figure,  the  oscillations  appear 
at  16.6  seconds  and  remain  until  the  end  of  the  flight.  Up  to  16.6  seconds, 
however,  the  values  of  Xn  are  well  determined.  In  order  to  obtain  mean¬ 
ingful  results  from  any  calculations  involving  Xw  from  16. 6  seconds  to  the 
end  of  the  flight,  the  oscillations  had  to  be  removed.  Since  the  scatter  shown 
in  Fig.  19  is  small,  it  was  felt  that  the  general  trend  of  \H  was  truly  re¬ 
presented.  Hence  a  curve  was  smoothed  through  the  fitted  values  as  shown 
in  Fig.  21.  The  results  of  the  fitted  and  smoothed  values  of  ‘Xm  >  which 
were  used  in  the  determination  of  the  aerodynamic  stability  derivative 
coefficients,  are  presented  in  Fig.  22. 

As  in  die  case  of  the  nutational  mode ,  the  Precession  Vector  dynamic 
damping  factors  were  well  determined  up  to  16.6  seconds  as  shown  in  Figs. 

23  and  24.  After  th.s  time,  however,  the  precessional  mode  began  to  be¬ 
come  extremely  .  majl  and  eventually  become  indeterminate  at  24.5  seconds 
from  release.  As  |Kp|  decreased,  was  more  difficult  to  determine 
accurately.  Phe  larger  scatter  of  Xp  shown  in  Fig.  23  and  the  increasing 
probable  error  shown  in  Fig.  24  indicate  that  the  results  were  not  representa¬ 
tive  of  the  general  trend  of  the  precessional  mode.  Thus,  to  smooth  the 
results  did  not  appear  reliable.  Therefore,  in  order  to  obtain  a  character¬ 
istic  dynamic  damping  factor  for  the  precessional  mode  from  16. 6  seconds 
tc  24.5  seconds  from  release,  the  magnitude  of  the  Precession  Vector  as  a 
function  of  time  from  release  was  analyzed.  ifffi  ,  too,  was  character¬ 
ized  by  oscillations.  Thus ,in  order  to  obtain  bona  fide  dynamic  damping  , 
factors,  as  a  function  of  time  from  release  was  smoothed  as  shown 

in  Fig.  25.  From  this  smoothed  curve,  the  natural  logarithm  of  IKpl  as 
a  function  of  time  from  release  was  obtained  as  presented  in  Fig.  26.  By 
taking  the  local  slopes  of  this  curve,  the  Precession  Vector  dynamic  damping 
factor  was  obtained  as  a  function  of  time  from  release  from  16.6  seconds  to 
24.5  seconds.  An  explanation  of  this  procedure  is  given  in  Appendix  I.  The 
combined  results  of  Xp  a  function  of  time  from  release  are  presented 
in  Fig.  27. 

It  should  be  noted  that  the  interpretation  of  the  results  of  die  fitting 
procedure  in  the  resonance  region  are  felt  to  be  questionable.  At  resonance, 
both  the  Nutation  and  Trim  arms  are  fixed  (in  body  axes),  and  the  ability  of 
the  fitting  technique  to  distinguish  between  the  two  is  uncertain.  For  the 
purpose  of  this  analysis,  however,  the  resulting  stability  parameters  were 
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interpreted  as  indicated  by  "Wobble”. 


Aerodynamic  Stability  Derivatives 

As  indicated  in  the  Computational  Procedure,  known  values  for  the 
normal  force  coefficient  derivative ,  C#*  »  are  necessary  for  calculations 

yielding  values  for  the  pitch  damping  moment  coefficient  derivative, 

Cm**Cm£  »  and  the  Magnus  moment  coefficient  derivative, Cm**  • 

Figs.  28  and  29  present  Cb*  as  a  function  of  time  from  release  and 
Mach  number,respectively.  The  wind  tunnel  test  results^  from  which  Cfe 
was  obtained,  are  presented  in  Appendix  I. 

In  order  to  determine  Cm*  as  a  function  of  time  from  release,  Eq. 
(22)  or  Eq.  (25)  may  be  used,  depending  on  whether  or  not  the  precessional 
mode  is  present.  One  method  is  to  use  Eq.  (22)  from  12.05  seconds  to 
24.5  seconds,  at  which  time  the  precessional  mode  is  no  longer  able  to  be 
determined  by  the  fit,  and  to  use  Eq.  (25)  from  24.5  seconds  to  the  end  of 
the  flight.  This  procedure  has  one  inherent  disadvantage,  however,  in  that 
for  a  few  seconds  beyond  24.5  seconds  the  precessional  mode  is  still  pre¬ 
sent  but  is  too  small  to  be  extracted  from  the  data;  thus,  a  slight  discontinuity 
will  appear  in  Cm*  as  a  function  of  time  from  release  at  24.5  seconds. 
Therefore,  in  order  to  avoid  this,  Eq.  (25)  was  used  to  determine  Cm* 
throughout  the  entire  flight. 

Using  values  of  ta**  as  a  function  of  time  from  release,  determined 
from  the  fitting  procedure,  and  the  roll  rate,  as  taken  from  Reference  4, 

Cm*  was  determined  as  a  function  of  time  from  release  by  Eq.  (25). 

These  results  are  presented  in  Fig.  30. 

From  these  values  of  Cm*  ,  the  gyroscopic  stability  factor,  s  ,  was 
determined  as  a  function  of  time  from  release  by  Eq.  (13).  These  results, 
presented  in  Fig.  31,  indicate  gyroscopic  stability  throughout  the  entire  flight. 
Near  resonance,  however, |ft  becomes  very  small.  This  is  a  direct  result 
of  the  large  values  of  Cm*  and  the  small  values  of  rolling  velocity  at  this 
time. 

An  examination  of  the  dynamic  damping  factors,  Figs.  22  and  27,  and 
die  gyroscopic  stability  factor ,  Fig.  31,  leads  to  the  result  that,  according 
to  the  linear  theory,  the  bomb  is  dynamically  stable  from  16.6  seconds  to 
25  seconds  from  release.  That  is,  it  is  only  during  this  time  interval  that 
both  necessary  and  sufficient  conditions  for  die  dynamic  stability  occur: 
X*vp<  0,  and  s  <  0. 

In  order  to  calculate  the  pitch  damping  moment  coefficient  derivative, 
Cmw+Cm*  »  b0**1  dynamic  damping  factors  must  be  known.  As  noted  pre¬ 
viously,  the  precessional  mode  was  not  able  to  be  extracted  from  the  data 
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after  24.5  seconds,  thus,Cww*CMi  may  be  calculated  from  Eq.  (23)  only 
up  to  24.5  seconds  from  release,  in  making  these  calculations  it  was  found 
that  the  pitch  damping  moment  coefficient  derivative  assumed  an  essentially 
constant  value  of  -49. 3  rad-*  from  20. 75  seconds  to  24.5  seconds  from  re¬ 
lease.  During  the  time  previous  to  this,  however, Cm®4  ranged  from 

-339. 3  rad"1  to  488. 3  rad  ,  changing  sign  at  approximately  14.2  seconds 
from  release.  Results  of  this  nature  are  unorthodox  to  the  linear  theory.  An 
examination  of  the  dynamic  damping  factors,  from  which  is 

calculated,  in  this  tint  period  indicates  that  these  results  would  indeed  occur 
since  both  \H  and  X?  change  sign  in  the  neighborhood  of  14  seconds  from 
release.  Furthermore,  calculations  of  the  Magnus  moment  coefficient  deriv¬ 
ative  from  12.05  seconds  to  24.5  seconds  from  release  from  Eq.  (24)  yielded 
results  of  a  similar  nature.  That  is,  previous  to  20.75  seconds  from  release, 
Cm*«  assumes  values  ranging  from  -642.2  rad"1  to  415.9  rad"1, 
changing  sign  at  approximately  14.2  seconds  from  release.  Both  of  these  re¬ 
sults  led  to  an  investigation  of  the  presence  of  some  other  phenomenon  at 
this  time  aside  from  resonance.  Again,  it  should  be  noted  that  these  results 
may  have  occurred  due  to  a  confusion  of  the  Nutation  and  Trim  arms  at 
resonance . 

In  the  Preliminary  Analysis,  it  was  noted  that,  from  12  seconds  to  20 
seconds  from  release,  the  rolling  velocity  assumed  a  fairly  constant  value, 
which  led  to  the  thought  that  Roll  Lock-In  and  Catastrophic  Yaw  effects  might 
also  be  present.  In  order  to  investigate  these  phenomena,  however,  a  know¬ 
ledge  of  the  roll  orientation  as  a  function  of  time  from  release  was  necessary. 

The  required  information  was  found  in  Reference  4.  Fig.  32  shows  the 
roll  orientation  angle,  ¥  ,  as  a  function  of  time  from  release.  It  was 
immediately  noted  that,from  12  seconds  to  approximately  21  seconds  from 
release,  the  bomb  experienced  the  effects  of  roll  orientation.  From  21 
seconds  to  the  end  of  the  flight,  however,  the  effect  was  cancelled  out.  For 
purposes  of  comparison,  the  magnitude  of  the  complex  angle  of  attack  as  a 
function  of  time  from  release  was  smoothed  and  is  presented  in  Fig.  33. 
Comparison  of  Figs.  32a  and  33a  shows  the  undesirable  effect  of  y  on 
H*l  up  to  approximately  21  seconds  from  release.  Correlation  of  the  roll 
orientation  angle  as  a  function  of  time,  Fig.  32a,  with  the  rolling  velocity  as 
a  function  of  time,  Fig.  6,  shows  that  the  misssile  was  on  the  threshold  of 
"lock-in"  in  the  vicinity  of  14  seconds  from  release.  Had  the  roll  moment 
due  to  cant  been  smaller,  or  the  induced  roll  moment  been  slightly  larger  in 
magnitude,  it  is  felt  that  the  bomb  would  have  experienced  Roll  Lock-In  and 
Catastrophic  Yaw.  A  small  change  in  the  flexible  tail  configuration  could 
well  induce  these  effects. 

Thus,  although  the  bomb  did  not  experience  Roil  Lock-In  and 
Catastrophic  Yaw  per  ae,  the  effects  of  the  induced  side  moment  associated 
with  these  phenomena  were  felt  up  to  approximately  21  seconds  from  release. 
Therefore,  the  unorthodox  results  obtained  by  applying  the  damping  equation 
of  the  linear  theory,  Eq.  (10),  to  the  data  from  12.05  seconds  to  20.75 
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seconds  from  release  were  attributed  to  the  presence  of  the  induced  side 
moment  during  this  time.  Hence,  the  dynamic  damping  factors  obtained 
from  the  fitting  procedure  from  12.05  seconds  to  20.75  seconds  from  re¬ 
lease  are  actually  X^jp  of  the  Catastrophic  Yaw  Theory. 

In  order  to  apply  the  resulting  equation  of  the  Catastrophic  Yaw  Theory, 
Eq.  (20),  certain  assumptions  had  to  be  made.  Fig.  34  shows  Cm»+Cm^  as 
a  function  of  time  from  release  calculated  from  the  linear  theory  from  16.7 
seconds  to  24.5  seconds  and  the  extrapolated  values.  As  seen  from  this 
figure,  CM«+CMi  was  assumed  constant  from  20.75  seconds  back  to 
12.05  seconds.  Xlso,  since  Xp  was  not  extracted  from  the  fitting  pro¬ 
cedure  beyond  24.5  seconds  from  release,  was  assumed  con¬ 

stant  from  this  point  to  the  end  of  the  flight.  The  linear  theory,  Eq.  (23), 
yielded  an  essentially  constant  value  of  Ciwu+Cm*  from  20.75  seconds  to 
24.5  seconds  from  release.  Thus,  in  extrapolating  Cm^+Cm*  on  both 
sides  of  this  time  interval,  the  assumed  value  was  found  to  be  the  same  on 
each  side,  -49.  3  .  It  is  lelt  that  this  assumption  of  constant 

was  justified  for  two  reasons:  1)  the  constant  value  obtained  from  the  fitting 
procedure  from  20.75  seconds  to  24.5  seconds,  and  2)  the  wind  tunnel  re¬ 
sults^  show  that,  although  the  moment  reference  position  is  different  than 
the  free  flight  case ,  +Cm*  varies  very  slightly  in  the  Mach  number 
range  experienced  in  free  flight.  See  Appendix  III.  Therefore,  the  pitch 
damping  moment  coefficient  derivative  was  obtained  as  a  function  of  time 
from  release  for  the  entire  time  of  flight  and  is  presented  in  Fig.  35.  From 
the  extrapolated  values  of  CMa+CM*  from  24. 5  seconds  to  the  end  of  the 
flight,  the  Precession  Vector  dynamic  damping  factor  during  this  time  period 
was  .calculated  from  Eq.  (10),  and  is  shown  in  Fig.  36  as  a  function  of  time 
from  release.  It  is  felt  that  the  fact  that  these  calculations  yielded  a  curve 
which  is  continuous  with  the  fitted  values  of  Xp  at  24.5  seconds,  and  which 
follows  the  trend  of  the  fitted  values,  enhances  the  justifications  made  in 
extrapolating  the  values  of  C*4+Cm*  • 

Another  assumption  was  necessary  concerning  the  Magnus  moment 
coefficient  derivative,  Cmp*  in  order  to  investigate  the  resonant  phase  of 
the  flight.  Fig.  37  shows  the  values  of  Cm**  obtained  from  the  linear 
theory,  Eq.  (24),  for  the  portion  of  the  flight  from  which  X**  and  X p 
were  extracted.  These  values  are  seen  to  increase  sharply  at  20.75  seconds 
from  release.  Thus, in  order  to  facilitate  investigation  of  induced  side  moment 
effects,  Cm**  was  assumed  constant,  having  a  value  of  13.5  rad‘l  , 
from  12.05  seconds  to  20.75  seconds  from  release.  This  extrapolation  is 
shown  in  Fig.  37.  From  24.5  seconds  from  release  to  the  end  of  the  flight, 
Cmp*  was  calculated  from  Eq.  (26)  using  the  extrapolated  values  of 
Cma+Cm*  These  results  are  shown  in  Fig.  38.  Upon  examination  of 
Eq7?10)  it  wa  j  noted  that  a  positive  Cm**  ^as  a  stabilizing  effect  on  the 
precessional  mode.  Hence,  after  resonance,  the  already  small  magnitude 
of  T?*  was  coupled  with  the  increasing  dynamic  pressure  and  the  stabilizing 
to  cause  the  precessional  mode  to  damp  out  halfway  through  the 
flight.  On  the  other  hand,  the  nutatioaal  mode,  which  is  predominant  after 
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24.5  seconds,  is  destabilized  by  the  positive  Magnus  term.  It  should  be 
noted,  however,  that  a  positive  Cw*  does  not  necessarily  imply  an  un¬ 
damping  of  the  Nutation  Vector.  The  combination  of  parameters  in  the 
Magnus  term  of  Eq.  (10)  must  be  larger  in  magnitude  than  the  combination 
of  parameters  in  the  pitch  damping  term.  In  this  case,  this  was  indeed  the 
situation.  Hence,  the  cause  of  the  undamped  motion  observed  in  the  latter 
half  of  the  flight  may  be  attributed  directly  to  a  Magnus  instability. 


As  soon  as  the  values  of  4 Ch*  and  Cwip<c  were  determined 

in  the  period  from  12.05  seconds  to  20. 75  seconds  from  release,  further 
investigation  of  the  Catastrophic  Yaw  effects  was  made.  Of  primary  interest 
are  the  dynamic  damping  factors  during  this  time  period,  which  are  pre¬ 
sented  in  Fig.  39.  It  should  be  noted  that,  since  Cm^  *  Cm*  is  directly 
proportional  to  the  mean  of  \H  and  \p  ,  and  inversely  proportional  to 
the  dynamic  pressure,  the  constant  value  of  Cm^Cm*  was  maintained 
as  a  result  of  the  increasing  value  of  Q.  The  fact  that  tne  observed  Magnus 
instability  was  slight  may  also  be  attributed  to  the  increasing  value  of  Q. 
The  values  of  the  dynamic  damping  factors  previous  to  20.75  seconds  are, 
adopting  Nicolaides'  notation,  p  .  Using  the  values  of  Cm^+Cm*  and 
Cm,*  between  12.05  seconds  and  20.75  seconds  from  release,  the  linear 
dynamic  damping  factors  may  be  calculated  from  Eq.  (10).  These  results 
are  shown  as  a  function  of  time  in  Fig.  40.  Fig.  41  presents  a  complete 
picture  of  the  effect  of  the  induced  side  moment  of  the  Catastrophic  Yaw 
Theory  on  the  dynamic  damping  factors. 


Employing  the  results  of  Fig.  41  and  Eqs.  (28a)  and  (28b),  where 
corresponds  to  the  linear  dynamic  damping  factors ,  yielded  the  induced  side 
moment  term  of  Eq.  (20)  as  a  function  of  time  from  release  as  shown  in  Fig. 
42.  From  Eq.  (19b)  and  the  pertinent  trajectory  parameters,  the  induced 
side  moment  derivative  was  calculated  and  is  presented  in  Fig.  43.  Eqs. 
(29a)  and  (29b)  yielded  CMy  ff,SIN4y  ,  which  is  shown  in  Fig.  44  as  a 
function  of  time  from  re  leaser 


It  is  immediately  obvious  from  the  results  of  Figs.  42,  43,  and  44  that 
the  calculations  from  the  nutational  and  precessional  modes  do  not  yield  the 
same  results.  This  is  a  direct  result  of  the  asymmetry  of  the  dynamic 
damping  factors  as  seen  in  Fig.  41.  Of  the  two  results  obtained,  it  is  felt 
that  the  calculations  involving  the  nutational  mode  are  more  indicative  of  the 
actual  phenomena  for  two  reasons:  1)  the  accuracy  to  which  the  Nutation 
Vector  dynamic  damping  factor  was  determined  was  significantly  greater 
than  that  of  the  precessional  mode,  and  2)  the  cause  of  the  termination  of  the 
fitting  routine  at  12.05  seconds  due  to  the  inability  of  the  iteration  pro¬ 
cedure  to  converge  to  a  value  of  |Kp|  and  .  As  can  be  seen  from  Fig. 
41,  the  Precession  Vector  dynamic  damping  factor  starts  to  decrease  at 
13.6  seconds  but  reverses  direction  at  13. 35  seconds  and  becomes  indeter¬ 
minate  at  12.05  seconds.  An  average  value  of  was  calculated 

from  values  obtained  from  the  nutational  and  precessional  modes.  For  the 
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reasons  mentioned  above,  the  average  value  is  not  felt  to  yield  better  results 
than  the  value  obtained  from  the  nutational  mode.  For  completeness,  how¬ 
ever,  the  average  value  is  presented  in  Fig.  45  as  a  function  of  time  from 
release. 

In  order  to  investigate  the  induced  side  moment  coefficient  as  a  func¬ 
tion  of  roll  orientation,  it  was  necessary  to  smooth  out  the  oscillations  in 
the  roll  orientation  as  a  function  of  time  from  release  curve,  Fig.  32.  The 
smoothed  curve  is  shown  in  Fig.  45.  Since  ¥  >£***._.  $t**4Y*  ,  and 
My  911441*  are  known  as  a  function  of  time  from  release,  cross -plotting 
yielded  Cm*.siN4V  and  My  sm4V  as  a  function  of  V  as  shown  in 
Figs.  47  ana48.  Again  the  most  reliable  result  is  felt  to  be  that  obtained 
from  the  nutational  mode.  It  was  noted  that  there  is  an  obvious  phase 
difference  between  the  resulting  curves  and  the  theoretical  curves.  The  theory 
predicts  that  the  period  be  W2  and  thatCM*u^H4V  and  have 

zero  magnitude  at  180°  and  270°.  Comparisonwith  wind  tunnel  results*®, 
howrrwr,  shows  that  the  period  observed  was  also  larger  than  *W/2 . 


Non-Linearities  in  Stability  Derivatives 

For  Round  739,  Reynolds  number,  Mach  number,  and  angle  of  attack 
were  varying  simultaneously.  Due  to  non-linear  effects,  it  was  not  possible, 
however,  to  determine  die  exact  functional  relationships  of  the  stability 
derivative  to  these  parameters  without  an  analysis  of  numerous  free  flight 
rounds. 

Angle  of  attack  variations  are  large,  from  7°  to  28°,  in  the  early 
portion  of  die  flight  while  small  variations,  from  7°  to  12°,  occur  during 
the  last  40  seconds.  Reynolds  number  also  varied  considerably,  ranging 
from  13x10®  at  the  beginning  of  the  flight  to  62x10®  at  the  end.  Mach  number, 
cn  the  other  hand,  experienced  a  relatively  small  variation,  from  0.75  to  0.95. 
It  should  be  noted,  however,  that  a  large  portion  of  the  flight  was  in  critical 
range  of  Mach  numbers,  from  0.90  to  0.95,Since  none  of  the  free  flight 
parameters  may  be  assumed  to  be  essentially  constant,  the  difficulty  of 
determination  of  relationships  was  not  decreased. 

Fig.  49  presents  Cm«  118  8  function  of  angle  of  attack  and  Mach 
number.  Substantial  increases  in  magnitude  can  be  seen  to  occur  with  in¬ 
creasing  angle  of  attack .  Also, Cm-  appears  to  be  non-linear  with  respect 
to  Mach  number,  even  though  the  change  in  Mach  number  is  smail.  it  should 
be  noted  that  since  die  fitting  procedure  obtains  the  stability  parameter  60^- 
from  which  Cm.  is  calculated,  by  fitting  over  a  prescribed  number  of  * 
cycles  of  data ,  the  angles  of  attack  used  for  this  analysis  are  obtained  by 
averaging  the  angle  of  attack  over  the  number  of  cycles  fitted.  Referring  to 
Figs.  30  and  33  it  can  be  seen  why  the  variations  were  observed  in  Cm« 
as  a  function  of  time.  increased  when  angle  of  attack  was  increasing 
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and  vice  versa.  An  exact  correlation  between  Figs.  30  and  33  cannot  be 
made  on  the  basis  alone,  however,  due  to  the  non-linearity  of  Cm,  with 
Mach  number.  A  curve  presenting  Cm*  as  a  function  of  Mach  number 
alone  would  not  be  representative  due  to  the  large  variations  of  angle  of 
attack  over  a  small  Mach  number  range.  It  should  also  be  noted  in  Fig.  49 
that  Reynolds  number  effects  are  also  present.  Although  one  would  expect 
the  restoring  moment  to  increase  as  Reynolds  number  increases  due  to  the 
delayed  separation,  it  is  felt  that,  coupled  with  the  fact  that  the  boat-tail 
effect  is  slight  and  the  magnitude  of  the  Reynolds  number  is  large ,  the 
variations  in  angle  of  attack  and  critical  Mach  number  are  more  influential; 
thus,  in  the  earlier  phase  of  the  flight,  Cm,  decreases  with  increasing 
Reynolds  number. 

Comparison  of  the  free  flight  results  with  those  from  wind  tunnel 
tests1^  indicate  the  free  flight  results  to  be  approximately  one-half  of  the 
wind  tunnel  values. Table  III  shows  typical  comparative  results.  The  wind 
tunnel  results  are  presented  in  Appendix  II.  Reference  3  indicates  that,  when 
analyzing  wind  tunnel  flow  visualization  results,  asymmetric  vortex  pairs 
were  observed.  It  was  further  mentioned  that  this  type  of  vortex  shedding 
is  subject  to  scale  effects  in  free  flight.  This  scaling  effect  presents  a 
possible  explanation  of  the  observed  difference  between  free  flight  and 
experimental  results. 

It  was  rot  possible  to  compare  Cm*  at  angles  of  attack  higher  than 
12°  due  to  the  fact  that  Cm  was  non-linear  with  oc  at  the  higher  angles 
of  attack^  and  thus  Cm*  could  not  be  computed  within  the  realm  of  the 
linear  theory.  It  was  possible,  however,  to  compare  Cm  •  This  is 
feasible  within  the  limits  of  the  linear  theory  due  to  the  fact  that,  in  the  fitting 
technique,  it  is  assumed  that  the  stability  parameters,  including  Cm*  » 
are  constant  over  each  sectional  fit.  Hence,  even  though  the  angle  of  attack 
during  the  sectional  fit  may  be  large,  Cm  may  be  obtained,  theoretically, 
by  multiplying  Cm*  by  the  corresponding  angle  of  attack,  assuming  zero 
trim.  Fig.  50  shows  these  results  obtained  compared  to  the  wind  tunnel 
results.  It  should  be  noted  that, in  comparisons  of  both  Cm*  and  Cm  , 
the  difference  between  free  flight  results  and  wind  tunnel  results  increases 
with  increasing  Mach  number.  It  is  felt  that  the  differences  observed  are 
due  to  possible  effects  of  the  critical  Mach  number  region.  In  Fig.  50,  at 
the  lower  Mach  numbers  the  difference  in  Cm  is  as  large  as  the  higher 
Mach  numbers,  but  this  is  primarily  due  to  the  fact  that  these  lower  Mach 
numbers  are  occurring  at  the  higher  angles  of  attack,  where  the  linear  theory 
assumption  used  to  calculate  Cm  Is  not  very  good.  But  at  the  lower  angles 
of  attack,  where  this  assumption  is  good,  one  would  expea  the  difference  to 
decrease.  It  does  not,  however,  and  it  is  felt  that  this  is  due  to  the  Mach 
number  becoming  more  critical.  Finally,  it  should  be  noted  that,  although 
the  magnitudes  of  the  parameter  are  different,  the  general  trends  are  the 
same. 
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TABLE  III.  COMPARISON  OF  FREE  FLIGHT  CM  TO  WIND  TUNNEL  VALUES 


Mach 

Humber 

Angle  of 
Attack  (deg) 

(1/rad) 
Free  Flight 

(1/rad) 

Wind  Tunnel 

% 

Difference 

0.80 

11.4 

-2.5 

-4.52 

44.2 

0.90 

10.5 

-2.3 

-4.45 

47.7 

0.95 

9.4 

-1.0 

-2.66 

62.5 
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-Cm  pitching  moment  coefficient 
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A  comparison  of  the  pitching  moment  coefficient  with  other  free  flight 
data  was  also  available.  The  results  of  Reference  18  agreed  very  well  with 
the  results  of  this  analysis.  While  the  Unified  Linear  Theory  does  not  include 
effects  of  roll  orientation,  it  was  expected  that  the  results  obtained  would  be 
averaged  over  roll  orientation.  Fig.  51  shows  the  free  flight  results  of 
Reference  18  and  of  this  analysis,  and  the  experimental  results  of  Reference 
16.  It  is  seen  that  the  results  of  the  Unified  Linear  Theoy  are  approximately 
the  average  over  the  roll  orientation  of  the  results  of  Reference  18. 

As  discussed  previously,  the  pitch  damping  moment  coefficient  deriva¬ 
tive  was  able  to  be  calculated  for  only  a  small  segment  of  the  flight  due  to 
the  induced  side  moment  effects  and  the  damping  out  of  the  Precession  Vector. 
The  values  calculated  were  essentially  constant,  and  thus  the  extrapolated 
values  were  also  equal  to  this  constant  value.  Comparison  with  wind  tunnel 
results  was  not  possible  due  to  the  fact  that  the  dynamic  tests ^  were  made 
with  a  center  of  gravity  position  different  from  the  free  flight  case.  The  free 
flight  c.g.  is  30.7%  of  the  body  length  from  the  ->ose  while  the  experimental 
c.g.  is  40%  of  the  body  length  from  the  nose.  Sii,:e  the  tests  were  made  at 
one  c.g.  position  only,  a  transformation  of  the  dynamic  moment  coefficient 
was  not  possible.  Although  comparison  of  the  magnitudes  was  not  valid,  a 
comparison  of  the  variation  of  C*u  as  a  function  of  Mach  number  in 
the  wind  tunnel  (see  Appendix  III)with  that  calculated  from  free  flight  data 
yielded  both  to  be  approximately  constant.  Thus,  as  mentioned  previously, 
the  assumption  of  extrapolating  C*w«<W  as  8  constant  equal  the  calculated 
amplitude  was  felt  to  be  a  valid  one? 

The  Magnus  moment  coefficient  derivative,  ,  yielded 

essentially  constant  results  throughout  the  entire  flight.  Fig. 5?  ‘presents 
88  8  fun0'*08  of  angle  of  attack  and  Mach  number.  As  discussed 
previously,  due  to  induced  side  moment  effects,  Cm**  was  assumed  con¬ 
stant  from  12.05  seconds  to  20.75  seconds  from  release.  As  Fig.  49  shows, 
this  time  interval  included  a  substantially  large  variation  in  angle  of  attack. 
Very  slight  variations  were  observed  with  Mach  number;  again,  it  is  noted 
that  the  critical  Mach  number  range, 0 .90  5  M  SO. 95,  is  characterized 
by  the  largest  variation.  The  general  continuity  of  the  extrapolated  values 
with  the  calculated  values  at  the  lower  Mach  numbers  ia  felt  to  add  assurance 
to  the  assumption  made  in  the  resonant  portion  of  the  flight.  In  accordance 
with  the  assumption  used  in  the  fitting  procedure, Cm*  may  be  calculated  as 
a  function  of  angle  of  attack,  as  presented  in  Fig.  53. 

Although  comparison  of  results  for  the  Induced  side  moment  coefficient 
derivative  was  difficult,  two  interesting  results  were  obtained.  Results  of 
wind  tunnel  tested  on  the  split-skit  bomb  yielded  that,  in  general,  the 
period  of  oscillation  of  the  induced  side  moment  coefficient  is  larger  than 
it t%  ,  as  predicted  theoretically.  Free  flight  results  indicate  this  same 

phenomena.  Although  free  flight  results  from  Reference  3  are  with  respect 
to  cruciform  finned  bombs,  the  body  profile  is  identical  to  the  split-skirt 
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bomb  (739).  These  results  indicate  that  the  peak  values  for  the  induced  side 
moment  vary  considerably  with  angle  of  attack,  and  in  general,  these  values 
occur  at  roll  orientation  angles  less  than  the  theoretical  values.  Although, 
in  the  free  flight  results,  constant  angle  of  attack  was  not  present,  the  peak 
value  of  the  induced  side  moment  coefficient  was  observed  to  occur  at  a 
smaller  value  of  roll  orientation  angle  than  the  theory  predicts. 
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SECTION  VI 


CONCLUSIONS 


Application  of  the  Method  of  Differential  Corrections  to  fit  the  Unified 
Linear  Aeroballistic  Theory  to  free  flight  angular  orientation  data  has  led  to 
a  determination  of  the  aerodynamic  moments  influencing  the  dynamic  stability 
of  a  split-skirt  bomb.  The  data  was  fitted  to  within  a  maximum  error  of  7% 
of  the  magnitude  of  the  motion.  The  aerodynamic  coefficient  derivatives 
which  were  determined  included:  pitching  moment  pitch  damping 

moment  (Cm^+Cm*  ),  ana  Magnus  moment  In  addition,  an  effort 

was  made  to  determine  the  induced  side  moment  coefficient  derivative, 

by  an  extension  of  the  linear  theory,  ^he  last  result,  however, 
was  questionable  due  to  complications  arising  in  the  fitting  technique.  The 
analysis  was  carried  out  assuming  no  such  complications  although  prelimi¬ 
nary  investigations  are  currently  being  carried  out  in  this  regard. 

The  stability  analysis  was  made  on  the  last  50  seconds  of  the  63-second 
flight.  During  this  time,  two  instabilities  were  determined:  Resonance  and 
Magnus .  Resonance  occurred  at  approximately  14  seconds  from  release 
while  Magnus  effects  were  felt  from  25  seconds  to  the  end  of  the  flight. 


Magnus  Instability 

At  25  seconds  from  release  the  precessions!  mode  damped  out  and  the 
remainder  of  the  flight  was  characterized  by  nutational  undamping  attributed 
to  a  positive  Magnus  moment.  During  this  phase  of  the  flight  the  motion  was 
essentially  circular,  increasing  in  amplitude  from  7°  to  approximately  12°. 

It  is  important  to  realize  that  the  cause  of  this  instability  does  not  lie 
in  the  destabilizing  Magnus  moment  alone.  From  the  Unified  Linear  Aero- 
ballistic  Theory,  it  is  seen  that  the  nutational  mode  does  not  necessarily 
undamp  when  a  positive  Magnus  moment  coefficient  derivative  is  present. 

Tli  j  Magnus  term  of  the  damping  equation  must  be  larger  in  magnitude  than 
the  pitch  damping  term  in  order  that  be  positive.  Each  of  these  terms 
is  composed  of  a  combination  of  mass  parameters,  trajectory  parameters, 
and  aerodynamic  coefficient  derivatives. 


Resonance  Instability 

At  approximately  14  seconds  from  release,  the  bomb  was  at  resonance. 
The  magnitude  of  the  complex  angle  of  attack  at  this  time  was  the  largest 
experienced  in  die  flight,  28°.  An  analysis  of  the  rolling  velocity  and  the  roll 
orientation  angle  as  a  function  of  time  showed  that  the  bomb  was  trying  to 
lock-in  for  approximately  6  seconds  in  the  region  of  resonance.  It  was  felt 
that  this  time  period  was  sufficiently  long  to  allow  the  induced  side  moment 


to  have  a  considerable  effect  on  the  magnitude  of  the  complex  angle  of 
attack.  Although  the  fitting  procedure  indicated  the  presence  of  an  induced 
side  moment  at  resonance,  further  investigations  are  currently  being  made 
to  explore  possible  numerical  problems  at  resonance.  This  results  from  a 
confusion  of  the  Trim  and  Nutation  arms  at  resonance.  It  was  felt, however, 
that  due  to  the  length  of  the  lock-in  region,  effects  of  the  induced  side 
moment  surely  were  present.  Hence,  although  the  results  presented  may  be 
magnified,  the  feeling  is  that  the  side  moment  is  generally  represented. 

Comparison  of  the  pitching  moment  coefficient  was  available  with 
both  wind  tunnel  and  other  free  flight  results.  The  wind  tunnel  results 
differed  considerably  from  results  of  this  analysis  while  the  other  free  flight 
results  compared  favorably.  It  was  felt  that  a  possible  explanation  for  the 
difference  between  free  flight  and  wind  tunnel  results  may  be  attributed  to 
scaling  effects  of  asymmetrical  vortex  shedding. 

Thus,  in  conclusion  it  was  demonstrated  that  the  fitting  technique 
developed  at  Notre  Darned  accurately  represented  the  angular  oscillation 
data  of  a  full  scale  bomb  in  free  flight.  From  these  fits,  the  pertinent  aero¬ 
dynamic  stability  coefficients  affecting  the  dynamic  stability  of  the  bomb 
were  accurately  determined.  The  analysis  indicated  that  the  bomb  experienced 
both  Resonance  and  Magnus  Instability.  Thus,  it  was  shown  that,  although 
the  split-skirt  tail  configuration  affords  greater  flexibility  in  retardation 
control  and  presents  a  more  streamlined  profile  to  the  flow  field  than  does 
the  cruciform -finned  configuration,  it  does  not  insure  freedom  from 
instabilities  characteristic  to  the  latter. 


APPENDIX  I 


LOG  DECREMENT  TECHNIQUE  APPLIED  TO  PRECESSION  VECTOR 


The  magnitude  of  the  Precession  Vector  may  be  expressed  as  a  function 
of  time  as  follows : 

t  .  ApH^te) 

iKpi-lKploC  (i-i) 

where  \Kf\o  is  the  initial  magnitude  at  time  tQ.  Taking  the  natural 
logarithm  of  Eq.  (I-l)yields 


(1-2) 


(1-3) 


Thus,  application  of  Eq.(I-3)  at  points  along  l-nltM  as  a  function 
of  time  from  release,  reinitializing  at  each  point,  yields  Xp  as  a  function 
of  time  from  release. 
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APPENDIX  II 


NOL  STATIC  WIND  TUNNEL  TEST  RESULTS 

Static  wind  tunnel  test  results**  from  the  Naval  Ordnance  Laboratory 
(NOL),  White  Oak,  were  used  to  calculate  the  normal  force  coefficient 
derivative ,  Czot  ,  as  a  function  of  Mach  number. 

Reference  1  presents  the  normal  force  coefficient,  Cg  ,  as  a  function 
of  angle  of  attack  and  Mach  number  with  roll  orientation  as  a  parameter .  A 
least  squares  fit  of  the  data  up  to  an  angle  of  attack  of  12  degrees  yielded 
Cg*  as  a  function  of  Mach  number,  with  roll  orientation  as  a  parameter, 
within  an  accuracy  of  1.0%.  These  results  are  shown  in  Fig. II- 1 .  Fig.II-2 
presents  the  same  results  averaged  over  roll  orientation.  It  is  these  values 
which  were  used  in  the  determination  of  the  aerodynamic  coefficient 
derivatives . 

A  similar  procedure  was  used  to  obtain  the  pitching  moment  coefficient 
derivative,  Cm*  ,  as  a  function  of  Mach  number.  These  values  were 
determined  for  the  purpose  of  comparison  with  the  free  flight  results.  In 
the  wind  tunnel,  the  model  was  tested  with  its  moment  reference  at  the  mid¬ 
point.  Hence,  for  comparison,  a  transformation  of  the  coefficient  was 
necessary.  The  transformation  equation  used  is  as  follows: 

Cmc<=Cm0( +  ^0*^  (II- 1> 

where 


*  refers  to  die  new  c.g.  positions 
XQ  is  the  distance  between  the  c.g.  position 
(positive  direction  is  forward) 
d  is  maximum  body  diameter 

Fig.  II- 3  presents  Cjvu  as  a  function  of  Mach  number  and  moment 
reference  (center  of  gravity)  position.  On  a  larger  scale,  Fig.  II-4 
Cm*  referred  to  the  free  flight  center  of  gravity.  It  is  these  values  which 
were  used  for  comparison  with  free  flight  results. 
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APPENDIX  III 


NOL  DYNAMIC  WIND  TUNNEL  TEST  RESULTS 


Dynamic  wind  tunnel  tests ^  from  the  Naval  Ordnance  Laboratory 
(NOL),  White  Oak,  yielded  the  pitch  damping  moment  coefficient  derivative 
as  a  function  of  Mach  number  and  roll  orientation. 

The  model  was  tested  with  the  center  of  gravity  located  40%  of  the 
body  length  from  the  nose.  Since  no  other  center  of  gravity  positions  were 
used,  a  transformation  of  c.g.  positions  was  not  possible.  For  completeness, 
however,  Cm  +  as  a  function  of  Mach  number  and  roll  orientation 
is  presented  ft  Fig.  Hi- 1,  while  Fig.  III-2  shows  Cmq  +  Cj^f  ,  averaged  over 
the  roll  orientation,  as  a  function  of  Mach  number.  "  * 
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13.  ABSTRACT  (continued) 


the  presence  of  an  induced  side  moment  at  resonance.  In  this  phase  of  the  analysis, 
however,  the  numerical  procedure  used  in  fitting  the  data  caused  concern  in  that, 
at  resonance,  neither  the  Nutation  Vector  nor  the  Trim  Vector  are  rotating  (i.e.  in 
body  axes).  Hence,  it  was  felt  that  in  fitting  this  portion  of  the  data,  the  numerical 
procedure  may  not  have  been  able  to  distinguish  between  the  two. 
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test  was  part  of  a  joint  research  proyram  on  instrumented  bombs  undertaken  by  the 
United  Kingdom,  Australia,  Canada,  and  the  United  States.  The  analysis  was  carried 
out  by  the  Department  of  Aero-Space  Engineering,  University  of  Notre  Dame.  The 
bomb  vas  observed  to  experience  Catastrophic  Yaw  in  the  early  phase  of  the  flight, 
while  a  Magnus  instability  was  evident  in  the  latter  phase.  An  analysis  of  the  last  hall| 
of  the  flight  dgta  yielded  an  excellent  determination  of  the  aerodynamic  stability 

C'am  a.  r.a  a  .  The  Unified  Linear  Aeroballistic  Theory  was  foum| 

Bll  during  the  Magnus  instability  phase  of  the  flight, 
ination  of  the  Magnus  stability  coefficient,  CM_. 
tigate  the  motion  of  the  bomb  as  it  passed  thrown 
served  at  this  time  coupled  with  a  substantial  in- 
>mplex  angle  of  attack,  thus  indicating  the  possibility 
tnd  Catastrophic  Yaw.  Fits  of  the  Unified  Linear 
Aerobaiiistic  Theory  to  the  data  Pere  successfully  accomplished  for  this  phase  of 
the  flight.  An  analysis  of  the  motion  yielded  stability  parameters  which  indicated 
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